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Abstract 

We propose a dual scale drift-diffusion model for interfacial growth and etch- 
ing processes. The two scales are: (i) a depletion layer width AW surrounding 
the aggregate and (ii) a drift length I. The interplay between these two an- 
tithetical scales yields a variety of distinct morphologies reported in electro- 
chemical deposition of metals, viscous fingering in fluids and in porous silicon 
formation. Further, our algorithm interpolates between existing growth mod- 
els (diffusion limited aggregation, ballistic deposition and Eden) for limiting 
values of AVFand 
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I. INTRODUCTION 



The study of complex patterns in dimension d = 2 has attracted a great deal of attention 
in the past decade and a half. Systematic studies of emergent patterns in electrochemical 
deposition (ECD) of metals and of the morphologies obtained during fluid-fluid displace- 
ment in a Hele-Shaw set up have been carried out On the other hand, complex pore 
geometries of anodically etched silicon have also evoked considerable interest, on account of 
the technological promise of porous silicon [Q. A simple model which would describe the 
fascinating varieties of morphologies obtained both in aggregation and dissolution would be 
a desirable objective. 

In electrochemical deposition experiments, two parameters, the electrolyte concentration 
(C) and the applied potential (V) are tuned to obtain (i) dendritic structures, both thick and 
needlelike (or stringy); (ii) Dense branching morphologies (DBM), which are homogeneous 
and (iii) randomly ramified self similar structures as in diffusion limited aggregation(DLA) 
0,^. In a Hele-Shaw cell, similar transitions are observed when the pressure and surface 
tension are varied In silicon, anodic etching gives rise to differing pore morphologies 
depending on the anodization potential and the substrate doping level 

A first principles explanation of such phenomena will have to encompass (i) the diffusive 
field, (ii) the Laplace field, (iii) convective processes and (iv) surface tension, curvature and 
underlying anisotropy effects. Hence, even numerical solutions may prove difficult. Simple 
formal approaches have been attempted for obtaining stability conditions at the interface 
It has also been hypothesized that the emergence of different characterizing length 
scales is the result of the interplay between the Laplace and diffusion fields governing the 
growth process in ECD [Q. Similar models have been proposed to model the anodic etching 
of silicon in hydrofiouric acid. A phenomenological Schottky diode model for porous silicon 
formation invokes a depletion layer width at the interface with a substantial barrier 
lowering at the pore tips. 

On the other hand, a simple algorithmic approach has been adopted in the last decade 
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and a half to successfully capture the basic patterns mentioned above p|-[TI|. The most 
common amongst these is the diffusion limited aggregation (DLA) [|12[ which results in a 
self similar structure resembling zinc leaves grown in an electrolyte under certain limiting 
conditions. Generalizations of the DLA have been proposed by the introduction of additional 
parameters into the model. Surface tension effects have been modeled by the incorporation 
of a sticking probability at the aggregate surface This resulted in a transition from 

the low density fractal clusters to regular patterns. Recently, a multiwalker DLA model 
has been proposed [|14| wherein two parameters, namely the particle concentration and the 
width of a "migrational envelope" are tuned to obtain dense branching morphologies as well 
as DLA patterns. 

In this article, we propose a single walker, two parameter drift-diffusion model to sim- 
ulate growth (and dissolution) in a wide variety of systems. Section 2 states the algorithm 
employed in the growth process. Section 3 summarizes the various growth patterns obtained 
in the simulations with differing parameters and correlates them to the morphological classes 
reported in experimental literature. Section 4 constitutes a brief discussion. 



II. THE DRIFT-DIFFUSION MODEL 

It will be a desirable goal to arrive at a model which can reproduce the various morpholo- 
gies obtained from different growth phenomena. Modeling the transition from one class of 
patterns to another will also be an interesting objective. Such transitions can be represented 
as resulting from the interplay of two or more parameters ||l5| , p!6[| . Different morphologies 
result in the various limiting cases. However, to maintain the simplicity of such an approach, 
a minimal set of parameters should be employed. 

In the classical DLA model [Q, a particle starts its random walk at an infinitely large 
distance away from a seed or cluster. Its random diffusive motion is terminated the moment 
it comes into contact with the aggregate cluster. This simulation models diffusion in the low 
concentration limit and does not employ any control parameters. To control the diffusive 
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process, the introduction of other variables such as a diffusion length [pT|,p!8| or particle 



concentration |[T^ have been suggested in literature. 

It is conceivable that very close to the aggregate surface, the particle movement is no 
longer controlled by the macroscopic diffusion field, but by microscopic, localized phenom- 
ena. To model such processes, additional parameters will have to be introduced into the 
simulation. We attempt to develop an algorithm which models pattern formation as a result 
of the interplay of a macroscopic field and localized surface phenomena. 

The main control parameters in this model are: (i) a depletion layer width APi^ which 
controls the diffusion and (ii) a drift length / governing the field driven processes in the 
proximity of the aggregate. The design of the algorithm is outlined below: 

(i) The particles are launched, one at a time, from a randomly chosen lattice site beyond 
the depletion layer boundary at a distance AW^from the surface of the aggregate. This 
is illustrated in Fig. 1. We begin with a circular depletion zone of radius APi^ around a 
central seed. As the aggregate grows, the depletion zone boundary is modified to follow the 
aggregate contour at a distance AW as shown in the figure. 

(ii) The particles execute random walks as in ordinary off-lattice DLA algorithms pi]| , pn 



(iii) The moment a part of the aggregate surface is encountered within a radius / of the 
particle location, the random walk is terminated. The particle is then moved to the surface 
(in a sense, "field driven" ) and becomes part of the aggregate. 

The simulations were stopped when the aggregate had grown to a radius of around 200 
units, one unit being the diameter of a single particle. Approximately 10^ — 10^ particles 
had to be launched for various cases. 

The two control parameters employed in this simulation have also been suggested else- 
where in literature. To model porous silicon formation, a "finite diffusion length" similar 
to our Alphas been suggested |T^. The drift length I in our simulation is in a way similar 



to the width of the "migrational envelope" of Erlebacher et al. [p!4|,pT|. A depletion layer 



width AW and a drift diffusion length / have been used in conjunction with other parame- 



ters to successfully model porous silicon formation p2|. Our algorithm interpolates between 



existing growth models in various limiting cases. In the limit / = 1 and AW oo, the 
algorithm is identical to the DLA. AW = and / = 1 approximates the Eden limit. As I 
becomes larger than AW (the case depicted in fig. 1), a fraction of the generated particles 
are directly transported to the surface without undergoing diffusion. This is akin to ballistic 
deposition. When / 3> AW, the ballistic process dominates the diffusive process. 

III. RESULTS 

Fig. 2 depicts the variety of morphologies obtained with varying Aiyand /. For small 
I /AW , we obtain DLA like patterns (Fig. 2a). A simple mass-radius scaling calculation for 
this pattern yielded a fractal dimension of 1.65 ± 0.05. Stringy structures similar to those 
reported in electrochemical deposition at high voltages are seen in the large Al^ limit 
with I /AW ^ 1 (Fig. 2b). The mass-radius scaling exponent fluctuates around ~ 1.45 
for this pattern. The stringy morphology observed in some experiments is essentially one 
dimensional in character This is reproduced for very large AW{= I) > 25. As AW{= I) 
is increased, the scaling exponent approaches unity (see flgure 4). On the other hand, for 
small AVTand I /AW ^ 1, we obtained short dense branches exhibiting a relatively smooth 
front which remains nearly circular throughout the growth period (Fig.2c). This has been 
identified in experimental literature as homogeneous or tip splitting patterns [3,S|7ip3| . The 
mass radius scaling for these structures resulted in a fractal dimension ^ 2. For 1/ AW > 3, 
thick dendritic growths with side branching was obtained (Fig. 2d). For large AW{> 5), 
the patterns become too inhomogeneous to show any well defined scaling However, for 
small AW, the model approaches the Eden limit and the resultant patterns are somewhat 
homogeneous. 

The occurrence of dendritic growth in such a simulation is surprising, since it is well 
understood that anisotropy is required in the interfacial dynamics for side branching to 
occur |T|. The present algorithm admits an underlying anisotropy in the way a particular 
surface site is chosen, when more than one point on the aggregate surface falls within a radius 
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/ of the random walker. This selection was done in three ways, (i) The nearest site was 
chosen, (ii) A site was chosen at random, (iii) The site nearest to the radial line connecting 
the random walker and the central seed was chosen. The morphologies remained similar 
for the case / < AW. The three methods yielded different results in the limit / ^ AW . 
Dendritic growth was observed only in case (iii), which corresponds to the imposition of 
a preferred radial field. In the case of electrodeposition, this could be interpreted as the 
applied radial electrostatic field. 

Figures 3 and 4 show the variation in the mass-radius scaling exponent df with varying 
values of AW and I. The values obtained for various AW when I is kept constant (Z = 2) 
are plotted in Fig.3. df is seen to vary from 2.00 to 1.65 ± 0.05 as AW goes from very small 
(~ 1) to large (> 15) values. This illustrates a smooth transition between Eden type growth 
and DLA. As AW is kept equal to / and varied, df decreases steadily with increasing AW 
(= /) and approaches unity. This is depicted in Fig. 4 and represents the emergence of 
stringy structures as shown in Fig. 2. 

Fig. 5 codifies our observations into a "phase" |Q plot. The plot was constructed on the 



basis of approximately a hundred patterns, grown to aggregate diameters of around 400 units. 
The transitions from one phase to another being continuous, the boundaries are not rigidly 
defined. For example, in the / ^ AW case, the dendritic patterns observed at large AW 
become denser and more homogeneous as AW approaches zero. The phases depicted have 
been reported earlier in electrochemical deposition experiments and related processes 

i- 

Several studies of two dimensional growth on a one dimensional substrate have been 
reported in the past. Matsushita et al. have grown zinc "trees" on a linear carbon cathode 
p5| . Such processes have also been modeled using a DLA like algorithm The etching 



process in porous silicon may also be viewed similarly [0,^,^. This is achieved in our 
model by using a linear substrate which is represented by a side of the lattice. The deposition 
is no longer radial, but "quasi-one dimensional". 

Fig. 6 depicts the dependence of the mean aggregate density on the drift length /. 

6 



Deposits were grown on a linear substrate of length 300 units for a given AW . A cross-over 
behavior is discernible, with a minima at AW = I. This can be understood in terms of our 
phase plot (fig. 5). The minima corresponds to the stringy region, where AW ~ /. On 
either side of the stringy phase, there exists phases of higher density. A similar behavior is 
seen in porous silicon where transitions are observed from a uniform network of thin pores 
to shorter, wider pores on increasing the doping level in the substrate On the other 
hand, increasing the applied potential leads to the formation of thick and relatively linear 
pipe-like pores 0. 



IV. DISCUSSION 

The finite diffusion length model of Smith and co-workers is known to exhibit a cross- 
over from fractal to non-fractal clusters when the aggregate size exceeds the diffusion length 
I8| , |27[] . In our model, a similar behavior is discernible for large aggregate sizes when the 



simulation is carried out for deposition on a linear substrate in a two dimensional lattice (see 
Figure 7a). The density decreases with depth and finally becomes almost depth independent. 
This is indicative of a non-fractal behavior (ie, in the depth independent region). However, 
in the case of radial deposition onto a central seed, no such transition was observed for 
aggregate diameters from 100 to 1000 lattice units. The mass radius scaling exponent 
remained independent of the aggregate size. This is depicted in figure 7b which plots the 
aggregate mass against the cluster radius. The constant slope of the log-log plot is indicative 
of scale invariance over a large range of diameters (100-1000). 

One can explore the phenomenological relevance of the parameters, AVTand /. The 
significant role played by microscopic dynamics in modifying the solutions of macroscopic 
analyses has received growing appreciation in the recent past |I| . In the diffusion field versus 
drift field dynamics governing growth in this simulation, AW is a control parameter for the 
macroscopic diffusion field, whereas / models microscopic surface effects. An example of 
a phenomenological model employing similar arguments is the Beale model [5] for porous 
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silicon formation which hypothesizes the existence of a depletion layer whose width varies 

as 



AVT oc 



(1) 



no 

where Va is the anodization potential, Vbi is a constant built in voltage and no is the 
substrate doping level. On the other hand, the parameter / can be related to the barrier 
lowering (A0s) due to microscopic surface irregularities. This lowering, A0s \/^s-, Eg 
being the enhanced local electric field. For planar interfaces, Es oc 0s/Al^ where 0s is the 
overall barrier height. Due to surface irregularities, this height is enhanced locally as 

As / — i> AW , Eg is very large at the tips, leading to stringy patterns and the minima in 



Fig. 6. In a recent work [^, it has been shown that AW and / can be correlated to the 
depletion layer width and barrier lowering given in equations (1) and (2). 

The existence of competing processes in electrochemical deposition and other phenom- 
ena is well established. The phase diagrams reported by Sawada et al. and Grier et al. 
represents various morphological classes resulting from a variation of the concentration and 
applied voltage. Assigning a direct correlation between our parameters and the experimen- 
tal parameters based on a comparison between the phase diagrams of the simulation and 
experiment is difficult. In all probability, AW and / are functions of both concentration and 
applied voltage. Further work needs to be done to establish the exact functional correlations 
between AW and / and the physical processes involved in electrochemical deposition. 

In a sense, the two scales AVTand / are antithetical, one separating the particle from 
the aggregate, the other driving it towards the aggregate. Various growth models (DLA, 
ballistic deposition, Eden) can be obtained as limiting cases of this algorithm. We stress that 
the two scales do not contain any a priori bias towards any specific morphological structure. 
Nevertheless, a variety of distinct morphological structures observed in experimental growth 
and dissolution is obtained in our simulation. 
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FIGURES 

FIG. 1. The drift-diffusion model. A particle is released beyond the depletion layer boundary 

(the closed curve) which dynamically follows the aggregate (solid circle cluster) contour at a dis- 
tance AW. When the particle wanders to within a radius I (see dotted arrow) of the aggregate 
surface, it is "field driven" to become a part of the cluster. 

FIG. 2. Map of the simulated growth patterns with varying depletion width ^AW ' and drift 
length (a) Patterns resembling DLA {AW = 10, I = 2), (b) stringy structures {AW = 10, 
I = 10), (c) homogeneous, dense branching morphologies with a nearly circular growth front 
{AW =2, I = 2) and (d) thick dendritic growth with side branches {AW = 8 and Z = 24). 

FIG. 3. The fractal dimension df plotted against AW (in lattice units) with I = 2. For small 
AW, df ~ 2, which is indicative of compact clusters. In the large AW limit, df stabilizes at 
1.65 ±0.05. 

FIG. 4. The fractal dimension df plotted against AW (in lattice units) with I equal to AW. 
df is seen to approach unity in the large AW{= I) limit. 

FIG. 5. "Phase" plot of patterns observed in the simulations. The phases depicted transform 
continuously into one another. Hence, the boundaries must not be treated as rigidly defining a 
transition. Both AW and I are expressed in terms of lattice units. 

FIG. 6. Plot of the percentage aggregate density vs. I (in lattice units) obtained in deposition 
on to a linear substrate of length 300 units. The two curves are for AW = 6 and AW = 9. In 
each case, a minimum is observed at AW = I. 

FIG. 7. (a). Percentage aggregate density vs depth (in lattice units) plot for the deposition on 
to a linear substrate. Simulation done for a 500x500 lattice with AW = 5 and / = 2. A cross-over 
to a non-fractal behavior where the density is almost depth independent is seen, (b) Log-log plot 
of the mass vs cluster radius for radial deposition on to a central seed. The clusters are seen to be 
scale invariant over a large range of diameters (100-1000 lattice units). 
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